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ABSTRACT 

We present near-IR images, obtained with the Hubble Space Telescope (HST) and the WFC3/IR 
camera, of six passive and massive galaxies at redshift 1.3 < z < 2.4 (SSFR< 1CT 2 Gyr" 1 ; stellar 
mass M ~ 10 11 M Q ), selected from the Great Observatories Origins Deep Survey (GOODS). These 
images, which have a spatial resolution of ~ 1.5 kpc, provide the deepest view of the optical rest-frame 
morphology of such systems to date. We find that the light profile of these galaxies is regular and well 
described by a Sersic model with index typical of today's spheroids. Their size, however, is generally 
much smaller than today's early types of similar stellar mass, with four out of six galaxies having 
r e ~ 1 kpc or less, in quantitative agreement with previous similar measures made at rest-frame 
UV wavelengths. The images reach limiting surface brightness fi ^26.5 mag arcsec -2 in the F160W 
bandpass; yet, there is no evidence of a faint halo in the galaxies of our sample, even in their stacked 
image. We also find that these galaxies have very weak "morphological k-correction" between the 
rest-frame UV (from the ACS z-band), and the rest-frame optical (WFC3 i?-band): the Sersic index, 
physical size and overall morphology are independent or only mildly dependent on the wavelength, 
within the errors. 

Subject headings: cosmology: observations — galaxies: fundamental parameters — galaxies: evolution 



1. INTRODUCTION 

Most of the stars observed in today's early-type galax- 
ies have formed at high redshifts (z > 2, see Renzini 
2006), but the observations have not yet constrained how 
that stellar mass has assembled in the presently observed 
systems. 

Several groups have reported that at least in the 
range M > M ~ 10 11 M©, the stellar mass func- 
tion of these high-redshift "elliptical galaxies", ap- 
pears to be rapidly increasing from z ~ 3 to z ~ 
1, pointing to this epoch as one of major assembly 
for modern bright galaxies (e.g. Bundy et al. 2006; 
Fontana et al. 2006; Arnouts et al. 2007; Scar- 
lata et al. 2007; Ilbert et al. 2009). A number of mech- 
anisms have been proposed as drivers of this evolution, 
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from merging to gradual in-situ accretion (Naab et al. 
2009; van Dokkum et al. 2009; Hopkins et al. 2010). 

Elliptical galaxies at z > 1.5 are often ob- 
served to be smaller, by factors ~ 3 — 5x, than 
their local counterparts of similar stellar mass, and 
thus to have » 30 — 100 x higher stellar density 
(Daddi et al. 2005; Trujillo et al. 2006; Trujillo et al. 2007; 
Toft et al. 2007; Zirm et al. 2007; Longhetti et al. 2007; 
Cimatti et al. 2008, Van Dokkum et. al. 2008, 
Buitrago et al. 2008, Toft et al. 2009), a fact that still 
lacks an explanation in terms of an evolutionary mecha- 
nism. 

Some have suggested that current observations at high 
redshift could have missed low surface-brightness halos 
surrounding the galaxies, due to the relatively limited 
sensitivity combined with the (1 + z) A cosmological dim- 
ming (Hopkins et al. 2009b; Mancini et al. 2009). Oth- 
ers have proposed that surveys of the local universe (e.g. 
the SDSS, Stoughton et al. 2002) could have missed a 
relatively large fraction of very compact, yet massive 
early-type galaxies due to the limited ground-based an- 
gular resolution (1" corresponds to a physical scale of 
~ 1 kpc at z = 0.05). While recent works (Valentin- 
uzzi et al. 2009; Trujillo et al. 2009) reported the identifi- 
cation of dense and massive galaxies in the local universe, 
their spatial abundance appears much smaller than the 
z ~ 2 galaxies. 

Furthermore, while small high-redshift samples have 
been imaged with HST and NICMOS at rest-frame opti- 
cal wavelengths (Buitrago et al. 2008; Trujillo et al. 2007; 
Toft et al. 2007; Zirm et al. 2007; Longhetti et al. 2007; 
Van Dokkum et. al. 2008), most of the galaxies have 
only been observed in the rest-frame UV with ACS, 
and the dependence of the morphology of these galax- 
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ies on the wavelength has not been characterized 
(Cimatti et al. 2008; Daddi et al. 2005). 

In this letter we use the unprecedented sensitivity and 
angular resolution of WFC3 /IR images recently acquired 
in the HUDF to study the optical rest-frame morphology 
of a small but well defined sample of low SSFR galaxies 
at z^2. AB magnitudes are used throughout this paper, 
and, when needed, a A-CDM world model with Ho=70 
km s" 1 Mpc" 1 , tt M = 0.3 and ft A = 0.7. 

2. DATA REDUCTION, SAMPLE SELECTION AND 
MORPHOLOGICAL ANALYSIS 

The primary imaging data consist of the first epoch 
of WFC3/IR observations acquired by Illingworth et al. 
(2009, program ID=G011563) in the HUDF (Beck- 
with et al. 2006), and thus the field is fully embedded 
in the GOODS-South field (Giavalisco et al. 2004). The 
images have been acquired in the F105 (Y), F125W (J) 
and F160W (H) filters and covers an area roughly equal 
to the footprint of the WFC3/IR camera (2.1 arcmin 2 ). 
We have carried out our independent reduction of the raw 
data and after rejecting images affected by persistence in 
the J-band, our final stacks reach 1 — a surface bright- 
ness fluctuations of 27.2, 26.6 and 26.3 AB/arcsec 2 in the 
three bands, respectively. We have drizzled the WFC3 
images to 0.06" pixel scale and brought them into reg- 
istration with the HUDF ACS and GOODS ones. The 
size of the PSF of the stacks is 0.18" in the H-band. 
We have then created a multi-wavelength source cata- 
log (GUTFIT, Grand Unified TFIT Catalog) using the 
TFIT procedures by Laidler et al. (2007), which contains 
PSF-matched photometry in all the bands acquired by 
the GOODS multi-facility program, from the U to the 
8^m IRAC one. 

We have selected a sample of high redshift, massive 
and passive galaxies using both the BzK color selection 
by Daddi et al. (2004), as well as multi-band SED fit- 
ting to spectral population synthesis models. While the 
BzK technique provides a well characterized criterion to 
select relatively massive galaxies with low specific star- 
formation rate at 1.4 < z < 2.5 , the method natu- 
rally has finite completeness and contamination. Thus, 
given the overall limited size of the sample, we have aug- 
mented it with galaxies at z > 1.3 selected for being mas- 
sive and passive according the criteria M> 10 10 M / M@ 
and SSFR< 10~ 2 Gyr" 1 based on fitting their observed 
broad-band SEDs to the Chariot & Bruzual (2007) and 
Maraston (2005) models. We have used Salpeter IMF 
and lower and upper mass limits 0.1 and 100 Mq, respec- 
tively, obtaining essentially identical results with both 
libraries. We have adopted the Calzetti et al. (2000) ob- 
scuration law to account for the possible (modest) pres- 
ence of dust. Finally, all our sample galaxies have no 
detection at 24/^m down to a 1 — a limit of 5/iJy, con- 
sistent with predictions for passively evolving galaxies at 
z ~ 2 (Fontana et al. 2009). 

The final sample consists of six galaxies, four 
of which have spectroscopic redshifts, namely # 
22704 and # 23555 from Cimatti et al. (2008), # 
24279 from Daddi et al. (2005) and # 24626 from 
Vanzella et al. (2008), and two for which we derived accu- 
rate (Az/(1 + z) ~ 0.05) photometric redshifts from the 
GUTFIT multi-band photometry. Four out of six galax- 
ies satisfy the pBzK criterion. Table 1 Q] lists the sample, 



including the parameters of the stellar populations from 
the fitted procedure, while Figure Q] shows the observed 
(points) and best -fit (curves) SED together with the H- 
band images of the galaxies. 

We have used the GALFIT package (Peng et al. 2002) 
to fit the light profile of the galaxies in the z—, Y-, J— 
and H- bands to the Sersic model 



I(r) = I e exp < -b n 



l/n 



1 



(1) 



where I(r) is the surface brightness measured at dis- 
tance r, I e is the surface brightness at the effective ra- 
dius r e and b n is a parameter related to the Sersic in- 
dex n. For n=l and n—A the Sersic profile reduces re- 
spectively to an exponential or deVaucouleurs profile. 
Bulge dominated objects typically have high n values 
(e.g. n > 2) and disk dominated objects have n around 
unity. Cassata et al. (2005), Ravindranath et al. (2006) 
and Cimatti et al. (2008) showed that GALFIT yields 
unbiased estimates of the Sersic index and effective ra- 
dius. 

We have obtained the PSF in each passband for use 
with GALFIT by averaging six well-exposed, unsatu- 
rated stars in the UDF/WFC3 field. For all galaxies, 
we ran multiple fits with the sky either set to a fixed, 
pre-established value or left as a free parameter. We also 
experimented with the size of the fitting region around 
each galaxy, finding peak-to-peak variations of the Sersic 
indices and r e at the 20% level at the most. The fits in 
Table Q] use a 6 x 6 arcsec 2 fitting region and free sky. 

3. THE MORPHOLOGICAL PROPERTIES OF GALAXIES 
AT Z ~ 2. 

Figure [5] shows the residual images obtained by sub- 
tracting the GALFIT best-fit models (after convolution 
with the PSF) from the corresponding images, together 
with the observed and best-fit model light profiles (see 
Table Q}. All galaxies, except # 24646, are very com- 
pact and symmetric, in close resemblance to present-day 
spheroidal galaxies. Galaxy # 24626, the most extended 
of our sample, has isophotes that deviate from the el- 
liptical shape, and the position angle of the isophotes 
varies as a function of the distance along the semi-major 
axis. These irregularities can be the result of ongoing or 
recent interactions (see Kormendy et al. 2009 and refer- 
ences therein). 

Of the six galaxies, four have Sersic index >2 (# 19389, 

# 22704, # 23555 and # 24626). In the case of ob- 
ject # 19389 we had to add a central unresolved (PSF) 
component to the Sersic model to obtain a reasonable 
fit. This central stellar object contains less than 10% 
of the total light in the object, and may be indicative 
of an AGN. This galaxy, however, is not detected in 
the Chandra Deep Field South 2-Megasecond Catalog 
(Luo et al. 2008) nor in the radio VLA maps by Keller- 
mann et al. (2008) and Miller et al. (2008). Figure[2]also 
shows that the inner part of the light profile of galaxy 

# 23495 (r< 1") is barely resolved in all bands, suggest- 
ing another AGN. This appears confirmed by the fact 
that the galaxy is also detected in the Chandra image, 
and it has X-ray luminosity 3.8±0.24 x 10 43 erg/s and 
5.6±0.6 x I0 43 erg/s in the soft and hard band, respec- 
tively; it is not, however, detected in the VLA maps. 
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Fig. 1. — The H-band images of the sample galaxies, and the stacked image of all but galaxy # 23495. The size of each panel is 
~ 3.5 X 3.5arcsec 2 . The insets show the observed (data points) and best-fit (curve) SED for each galaxy. For object # 23495 we also show 
the AGN SED best-fit as a dashed curve (see text). 
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Fig. 2. — The H-band residual images of the sample galaxies and of the stacked image obtained subtracting the best-fit Sersic model 
(after PSF convolution) from the original images. The insets show the observed profiles (diamonds), the best-fit profiles (continuous line), 
the psf profile (dotted line), a Sersic model with r e = 3kpc and ng eraic = 4 (dashed line) and the residual profile (bottom panel). The sky 
level is shown as an horizontal line. Black symbols refer to the H-band light profile, red diamonds to the Y— band one. The red Y-H colors 
of our galaxies cause the Y-band profile to be always below the H-band one. 



With the exception of the IRAC 8 /im data point, whose 
value suggests the presence of a hot dust component, the 
broad-band SED of this source is otherwise very similar 
to that of the other galaxies, and fits made under the 
assumption that its luminosity is powered by stars yield 
"stellar populations parameters" that are essentially the 
same as the other galaxies (see Table [T]) . We also note 
that fitting the SED of this galaxy with the AGN tem- 
plate by Polletta et al. (2006) returns the same value of 



the photometric redshift as with stellar templates. 

Galaxy # 24279 has a Sersic index of intermedi- 
ate value n^l.6, suggesting that this is probably a 
bulge dominated disk galaxy. We note that Stock- 
ton et al. (2008) and McGrath et al. (2008) have found 
that some passive galaxies at z > 1.5 have have disk 
morphology. Finally, to obtain a reasonable fit to the 
light profile of galaxy 24626, we had to use a combi- 
nation of two components: one with a very high Sersic 
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TABLE 1 

The sample of passive galaxies 



ID 


RA 


DEC 


zW 


log(A//M ) 


log(SSFR) 


E(B-V) 


Age 


T 


Sersic 


r e 












[Gyr- 1 ] 






Gyr 


//—band 


[kpc] 


19389 


53.135730 


-27.784932 


1.307p 


10.41 


-5.55 


0.1 


3 


0.2 


2.99 ( 2 ) 


1.02 W 


22704 


53.153799 


-27.774587 


1.384s 


10.70 


-5.55 


0.15 


3 


0.2 


2.65 


0.50 


23495 


53.158455 


-27.773982 


2.349p 


11.14 


-3.39 


0.15 


2 


0.2 


(3) 


< 0.38 


23555 


53.158810 


-27.797155 


1.921s 


10.82 


-2.98 


0.0 


2 


0.2 


1.97 


0.44 


24279 


53.163005 


-27.797655 


1.980s 


10.63 


-3.39 


0.15 


3 


0.5 


1.63 


0.37 


24626 


53.165159 


-27.785869 


1.317s 


11.10 


-2.11 


0.1 


3 


0.5 


7.42 


3.69 



Note. — (1): s and p indicate spectroscopic and photometric redshift, respectively; (2) modeled with the Sersic profile plus the psf one; 
(3) modeled with the psf profile, so the size reported here is just an upper limit. 



index to reproduce the inner part of the galaxy, and a 
second with n^l to reproduce the outskirts. Further- 
more, both components have effective radius r e ~3 kpc, 
implying that even if the component with lower surface 
brightness were not detected (e.g. in shallower images), 
the effective radius would still be found to be ~ 3 kpc, 
significantly larger than the other galaxies. Cases like 
this galaxy, therefore, do not provide an explanation to 
the much smaller size of the z ~ 2 ellipticals compared 
to those at z ~ 0. 

We have also analyzed the residual maps produced by 
subtracting the best-fit models from the corresponding 
original images to investigate the presence of diffuse low- 
surface brightness halos, as recently suggested by some 
(Mancini et al. 2009; Hopkins et al. 2009b). For all galax- 
ies, these maps (see Fig. [5]) show the presence of residual 
structure in the innermost 0.5" of the light profile, where 
the contribution by the PSF is the largest. This is very 
likely the result of the variations of the PSF across the 
field, which cannot be taken into account by our average 
PSF. To verify this possibility we have fitted the individ- 
ual stars using the same average PSF used for galaxies, 
and indeed we have found residuals similar in intensity 
and morphology to those of the galaxies. There is, how- 
ever, no evidence of a halo. The only exception is the 
case of galaxy # 24626, which requires a two-component 
model to reproduce its light profile, an inner bulge and 
an outer disk. Subtracting both of these two compo- 
nents, however, leaves a complex residual map, revealing 
a substantially more irregular light distribution than the 
other galaxies. 

To further explore the possible presence of a halo, 
we have stacked all the individual images together, ex- 
cept galaxy # 24626. In the H-band the stacked image 
has an equivalent exposure time T exp = 78, 500 sec and 
reaches la surface brightness sensitivity ^h, stack ^27 
mag arcsec" 2 , and we have performed the GALFIT anal- 
ysis on it following the same procedure used for individ- 
ual galaxies. The stacked images and the GALFIT resid- 
uals are shown in Figures [T] and [2J respectively. Close 
companions were not masked before stacking and were 
not included in the GALFIT model. The stacked image, 
like the individual galaxies, is very compact. The resid- 
ual images for individual galaxies and the stack show no 
evidence of a large-scale diffuse halo with less than 2% 
of the residual light falling within 2" of the source. 

We plot in Figure [3] the Sersic index (top) and effec- 
tive radius (bottom) of the galaxies measured from the 
the z, Y, J and H-band images as a function of wave- 
length. Assuming a mean redshift < z >= 1.7 for the 
sample, the covered rest-frame wavelength range extends 




Fig. 3. — The Sersic index (upper panel) and the effective radius 
(bottom panel) as a function of the wavelength. Object # 23495, 
the one best-fitted by a point like source, is not shown in this figure. 
The typical error bar is marked with a vertical bar, and amounts 
to ~ 30% of the measured Sersic indices and ~ 20% of r e . 



from the rest-frame UV at A ~ 3300 A to the optical at 
A ~ 6000 A. Consistent with a visual inspection, the two 
plots show weak or no morphological k-correction. The 
Sersic index is constant with wavelength within the er- 
rors, while the effective radius of three galaxies has a 
weak, but statistically significant dependence on it, de- 
creasing by about ~ 40% from the z band to the H one 
(r e is constant for the reminder of the sample). This 
weak dependence of r e with wavelength results in a neg- 
ative color gradient, with the outer parts of the galax- 
ies being bluer than the center. This color gradient has 
also been observed by Guo et al. (in preparation), who 
report that SED fitting implies that the stellar popula- 
tions in the outer regions of these galaxies are, on aver- 
age, younger by ~ 0.5 Gyr than those in the center. Fi- 
nally, we note that the observed weak dependence of mor- 
phology on wavelength is in good quantitative agreement 
with other studies of early-type galaxies at high-redshift 
from ACS and NICMOS images (McCarthy et al. 2007; 
Trujillo et al. 2007) as well as at lower redshift (Pa- 
povich et al. 2003; Cassata et al. 2005). 

4. THE MASS-SIZE RELATION AT Z ~2 

In Figure [4] we compare the relationship between stel- 
lar mass and size in the H band for the six galaxies in 
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Fig. 4. — The effective radius as a function of stellar mass for 
early— type galaxies. Our sample is marked with red filled circles, 
with galaxy # 23495 shown as an upper limit, since it is unre- 
solved. Crosses, circles, clovers, squares, stars, triangles, pentagons 
and diamonds represent the galaxies from Cimatti et al. (2008), 
Zirm et al. (2007), Toft et al. (2007), Daddi et al. (2005), Van 
Dokkum et al. (2008), Trujillo et al. (2006), Mancini et al. (2009) 
and Longhetti et al. (2007), respectively. The continuous and 
dashed lines represent the local relation and its 1 — a scatter, from 
Shen et al. (2003). 



our sample to other measures from the literature, both at 
high redshift and in the local universe. When needed, we 
have homogenized mass measures to M05 and CB07 with 
Salpeter IMF by applying offsets^ computed by Salim- 
beni et al. (2009). The figure shows that four out of our 
six galaxies are located below the local relation for early- 
type galaxies. The remaining two galaxies (#19389 and 
#24626) are on the local relation. 

Thus, on the one hand our small sample of ~ 2 massive 
and passive galaxies seems to confirm that a significant 
fraction of these systems are ~ 3 to 5 times smaller, and 
hence have ~ 50-100 times denser central stellar den- 
sity, than their local counterparts, as previously reported 
from similar but shallower and/or lower quality images 
(Trujillo et al. 2007; Toft et al. 2007; Zirm et al. 2007; 
Longhetti et al. 2007; Van Dokkum et al. 2008; 
Toft ct al. 2009) or at UV rest-frame wavelengths 
(Daddi et al. 2005; Cimatti et al. 2008). On the other 
hand, it also suggests a diversity of morphological proper- 
ties (size and stellar mass density) among these galaxies, 
with some of them being similar to the local ones, which 
still remains to be characterized from a statistical point 
of view. We plan to report soon on one such study. 

5. SUMMARY AND CONCLUSIONS 

Massive galaxies (M ~ 10 11 M Q ) with early spectral 
type at z ~ 2 appear to often have a relatively smooth 
and regular morphology, nearly independent of the wave- 
length, and with light profile which is well described by 
a Sersic model, i.e. their morphology can be classified as 
"elliptical" , although systems with an exponential com- 
ponent are also observed. This suggests that a tight cor- 



11 log(M C S07,Af05) 

log(M chab ) +0.25 



log(M SO 03) - 0.2 and log(M Saip ) 



relation between spectral and morphological properties 
similar to what is observed in the local universe, i.e. the 
back bone of the Hubble sequence, is already in place 
among bright galaxies at ~ 25% of the cosmic time. The 
Sersic analysis has also quantitatively shown that the 
morphology of these galaxies depends very mildly on the 
wavelength, between the rest-frame UV and optical win- 
dows. 

One of our six galaxies is almost certainly an AGN, 
while another with a spatially unresolved nuclear com- 
ponent is a candidate one. Such kind of contamination 
should not be surprising given the similarity between the 
SED of obscured AGN and quiescent galaxies. By the 
same token, however, it is also possible that both types 
of sources contribute to the luminous output of the same 
galaxy, a fact that we have not directly investigated in 
this work. The apparently high level of contamination 
(33% of our sample) observed at z ~ 2, the cosmic era 
when both AGN and star-formation activity both reach 
their maximum (Ueda et al. 2003), raises the intriguing 
question of whether this is direct evidence of a link be- 
tween the AGN and the mechanisms responsible for the 
cessation of star formation in the galaxies. We plan to 
report on this issue using larger samples from imminent 
new observations with WFC3 of the GOODS fields. The 
high AGN contamination is also likely to bias statistical 
conclusions made from UV/Optically-selected samples 
of passive galaxies at z ~ 2, unless sensitive X-ray and 
mid-IR photometry (Daddi et al. 2005, 2007), or at least 
multi-wavelength high-resolution HST imaging, can be 
used (like in the GOODS) to cull these sources. In sur- 
veys where such information is not available, like some 
recent ground-based ones (e.g. van Dokkum et al. 2009), 
the AGN so frequently encountered among the galaxies 
cannot be recognized. 

Four out of six of our z ~ 2 ellipticals lie below the local 
mass-size relation for early-type galaxies, in agreement 
with previous findings that they can be up to ~ 3-5 times 
more compact at given mass, and thus ~ 30-100 times 
denser, than their local counterparts. Of the remaining 
two galaxies, however, one (# 19389), which has an unre- 
solved central component, is at the boundary of the local 
relation, while the other (# 24626) is well within it. We 
have not found any evidence of a low-surface brightness 
halo surrounding these galaxies, either in the individual 
images, or in their stacked one. At most, only < 2% of 
the light (and hence of the mass, assuming a constant 
M/L ratio) can be segregated in such a halo. Galaxy # 
24626, requires a two-component light profile, one Sersic 
model with high n (spheroid) embedded in a disk, to 
provide a good fit to the observations, but both have the 
same, relatively large, effective radius, r e ~ 3 kpc. Thus, 
what is emerging here is that these galaxies are charac- 
terized by a dispersion of properties that also includes 
systems similar to the local ones, as some have already 
suggested (Cimatti et al. 2008; Mancini et al. 2009), 
which however still needs to be studied and character- 
ized. Until this work is done, whether and/or how the 
ultra-compact z ~ 2 ellipticals evolve into the local ones 
or whether ultra-compact systems are also present in 
the local universe but have so far eluded detection will 
remain open questions. 
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